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ABSTRACT 



We report on K-band VLTl/ AMBER observations at medium spectral resolution (~1500) of RS Capricomi, an M6/M7111 semi-regular 
AGB star. From the spectrally dispersed visibilities, we measure the star diameter as a function of observing wavelength from 2. 13 to 
2.47 microns. We derive aRosseland angular diameter of 7.95 + 0.07 mas, which corresponds to an effective temperature of 3160 ± 160 
K. We detect size variations of around 10% in the CO band heads, which are indicative of strong opacity effects of CO in the stellar 
photosphere. We also detect a linear increase in the size as a function of wavelength, beginning at 2.29 microns. Models of the stellar 
atmosphere, based on the mass of the star estimated from stellar-evolution models, predict CO-size effects of about half the sizes 
observed, and cannot reproduce the linear size increase with wavelength, redward of 2.29 microns. We are able to model this linear 
size increase with the addition of an extended water-vapor envelope around the star. However, we are unable to fit the data in the 
CO bandheads. Either the mass of the star is overestimated by the stellar-evolution models and/or there is an additional extended CO 
envelope in the outer part of the atmosphere. In any case, neither the water-vapor envelope nor the CO envelope can be explained 
using the current models. 
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' The effective size of a star at a given wavelength depends on 
. the opacity of the stellar atmosphere at that wavelength. Since 
' we effectively measure the diameter of the t - 1 surface of the 
star, the size is related to the extension of the atmospheric region 
where the absorption is produced. The diameter of that surface 
varies with opacity, hence frequency. The atmospheres of cool 
giants are so extended that these size variations are observable 
with an interferometer Quirrenbach et al. (1993 2001) studied 
the extended absorption regions of stellar atmospheres of cool 
. giant stars in the TiO band at 712 nm with the Mklll interferom- 
■ eter. These authors measured the visibilities of a set of 47 stars 
" using two filters, one centered in the TiO band, and the other 
in the continuum part of the spectrum close to the TiO band 
(754 nm). They found that the sizes in the TiO band are larger 
, than those in the continuum, and that this effect is stronger for 
cooler stars. The sizes in the TiO absorption band are ~ 10% 
larger than in the continuum for R - I color indices of ~ 1.6 
(spectral types M3-M4) and as much as 30% larger for R - I 
color indices of ~ 2.2 (spectral types M6-M7). 

Similar differences between the sizes observed in the con- 
tinuum and in bands containing absorption band heads of other 
molecules (such as H2O or CO) have been reported for AGB 
stars. Mennesson et al. (I2002I I and Perrin et al. (I2004I I. for in- 
stance, reported very large size ratios (50% or more) for some 
stars. 



* Based on observations made with ESO Telescopes at the Paranal 
Observatory under programme ID 383.D-0619 



Quirrenbach et al. ( 119931 1200 11 1 were successful in quali- 
tatively reproducing their data with the latest set of cool gi- 
ant models from the general-purpose stellar atmosphere code 
PHOENIX. Spherical, hydrostatic, massively line-blanketed at- 
mosphere models were constructed and used to predict the 
uniform-disk diameters in the TiO band and the continuum band 
as a function of model effective temperature, surface gravity, and 
mass (the stellar mass was used in the modeling, since it con- 
trols the deviation from plane-parallel atmospheres.) For most 
of the observed oxygen-rich giants, the diameter ratios of the 
TiO band to the continuum band agreed with the models com- 
puted for masses ~ 0.5 Mq. Hence, model atmospheres with very 
low stellar masses could fit the large diameter ratios observed in 
many stars, although evolutionary models predict masses as high 
as 5 Mq. 

A possible explanation of this inconsistency would be the 
existence of a transition zone at the base of the stellar wind 
(the MOLsphere; e.g. Tsuji 2008 and references therein), which 
could provide sufficient opacity in the TiO bands (and other 
molecular bands) to make AGB stars appear much larger 
than predicted by the hydrostatic model atmospheres such as 
PHOENIX. According to this picture, one would expect similar 
size effects for the CO band heads. Therefore, we decided to use 
AMBER in the K band to measure the effective sizes of a set 
of four cool giant stars through the CO band heads at 2.3 pm. 
The use of AMBER in medium-resolution mode (A/AA ~ 1500) 
provides considerably more information than could be obtained 
with the Mklll interferometer (Quirrenbach et al. I1993I I and the 
IOTA interferometer (Mennesson et al. 120021 Perrin et al. l2004l l. 
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in which narrow-band filters were used. In this paper, we report 
on the resuhs obtained from the analysis of the observations of 
the first star of our sample: RS Cap. 

RS Cap (HD 200994) has a K-band magnitude of -0.2 
(Cutri et al. 120031). It is a semi-regular variable (SRb) of spec- 
tral type M6/M7III and is located at a = 2 Ih 07m 15.4s, 6 - 
-16° 25' 21.4" (J2000.0). It has a visual magnitude of 8.3 and a 
parallax n - 1 .26 + 0.86 mas (van Leeuwen 2007 ), which maps 
into a distance between 500 and 2500 pc. However, following 
Scalo (1 19761 1 or Winters et al. (I2003I I. a bolometric absolute mag- 
nitude between -4 and -5 is derived for RS Cap (around -4.8 
in the case of Winters et al.). These estimates, in addition to the 
spectrum fitted between 0. 1 and 2jum, translate into a distance of 
-300 pc (Richichi et al. [T992] l. 

A variability amplitude of AV ~ 0.5 is seen in the Hipparcos 
data (Ferryman & ESA 1997 ), although a larger variability am- 
plitude (in photometric magnitude) of AB ~ 2 was reported in 
Kukarkin et al. (1969) , which has a period of ~ 340 days. 

The remainder of this paper is structured as follows: In Sect. 
|2] we describe our AMBER observations and the strategy fol- 
lowed in the data calibration and reduction. In Sect. [3] we report 
on the results obtained: Sect. 13. H is centered on our estimate of 
the effective temperature and Sect. 13.21 on the size effects ob- 
served at the CO band heads. In Sect.|4] we compare our results 
with synthetic data obtained from model stellar atmospheres. In 
Sect.|5] we summarize our conclusions. 

2. Observations and data reduction 

We observed RS Cap with the ESO Very Large Telescope 
Interferometer (VLTI) using the Astronomical Multi-BEam 
combineR, AMBER (see Petrov et al. 2007 for details on this 
instrument), in medium-resolution mode (AjAA ~ 1500). This 
instrument performs simultaneous observations of the inter- 
ferometric fringes generated by three telescopes. Therefore, it 
measures closure phases, which are quantities independent of 
atmospheric or instrumental telescope-dependent contributions 
(e.g. Rogers et al. |1974| l. AMBER also measures the so-called 
differential phase on each baseline, which roughly represents 
the phase of spectral features with respect to that in the con- 
tinuunfl The observations were performed on 4 June 2009, 
from 06:30 UT to 08:30 UT, using the three VLTI Auxiliary 
Telescopes (AT) positioned on stations EO, GO, and HO. These 
stations are distributed roughly in the east-west direction. In 
Table [1] we give the projected baseline lengths and position an- 
gles for our observations, together with a summary of the at- 
mospheric observing conditions. 18 Cap (MOIII star located at 
a = 20h51m49.3s, S = -26° 55' 08.9", J2000.0) was also ob- 
served as a calibrator, with the same observing configuration 
used to obtain the RS Cap visibilities. 

The high flux densities of RS Cap and 18 Cap in the H band 
(magnitudes of 0.25 and 0.42, respectively), allowed us to use 
the fringe tracker FINITO (see Gai et al. I2004I I. which employs 
part of the H-band light to correct, in real time, the delay shifts 
of the fringes due to atmospheric turbulences, thus increasing 
the coherence of the signal. A larger coherence allows the use of 
a longer detector integration time (DIT). Given the moderately 
good atmospheric conditions during our observations (see Table 
[T]l, we achieved a DIT of 200 ms with a negligible signal loss due 

' To be precise, two pieces of information are lost in the differential 
phase, compared to the original Fourier phase: the phase offset (or phase 
at the first spectral channel) and the phase gradient (phase as a linear 
function of wavenumber). 



to atmospheric jitter This allowed us to simultaneously observe 
a large spectral range of the K-band (between 2.13 and 2.47 yum) 
using medium-resolution mode. This wavelength range contains 
the i^co (2-0), (3-1), and (4-2), as well as the '^CO (2-0) 
band heads, together with the continuum blueward of the '^CO 
(2-0) band head at 2.29 fim. 

The data acquisition was carried out in the following order: 
(1) observation of an artificial signal for calibration of the in- 
strumental dispersive effects; (2) one exposure of dark frames 
(an exposure consisted on 200 frames of 200 ms each); (3) five 
exposures of the target; and (4) one exposure of empty sky close 
to the target. The observations of dark, target, and sky exposures 
were iterated once for 18 Cap (spanning 10 minutes), three times 
for RSCap (spanning 1 hour), and one more time for 18 Cap 
(spanning 10 minutes). The dead time between iterations, when 
the sources were not observed, was dedicated to the setup of 
the instrument: preparation of the optical-delay lines, telescope 
pointing, etc. 

The raw visibilities and differential closure phases for RS 
Cap and 18 Cap were obtained using the amdlib libraries (ver- 
sion 2.2) and the interface provided by the Jean-Marie Mariotti 
Center (JMMC). First, we manually aligned the spectrally- 
dispersed photometry channels of the ATs to the interferomet- 
ric channel, based on observations of the target and calibrator 
stars. We removed the bad pixels and took into account the flat, 
dark, and sky contributions. Afterwards, we calibrated the in- 
strumental dispersive effects and fringe-fitted each frame of RS 
Cap and 18 Cap (based on the P2VM algorithm of Tatulli et al. 
120071) . The resulting visibilities of individual frames were se- 
lected and averaged for each exposure in several ways, to test 
the robustness of our results on different averaging and selec- 
tion schemes. The finally chosen selection scheme was based on 
an atmospheric piston criterion (keeping only frames with a pis- 
ton smaller than 8 //m, to select only well-detected and centered 
fringes in all three baselines) and keep 50% of the remaining 
frames based on a signal-to-noise ratio (SNR) criterion. 

After the frame selection, we obtained a single frame- 
averaged visibility spectrum for each exposure. The last step in 
the data reduction was to average the visibilities of all exposures 
of each star to increase the SNR. This step was performed out- 
side amdlib, with an in-house developed python-based program 
that uses the PyFITS Ubrary (provided by the Space Telescope 
Science Institute, operated by AURA for NASA). The averaged 
closure phases and differential phases (i.e., those resulting from 
the averaged exposures) are zero within 3°, in absolute value, for 
both stars and for all the spectral channels. These small values of 
the closure phases indicate that the emission from both sources 
is symmetric, at least in the direction of the projected baselines 
of our observations and to the scale corresponding to our spatial 
resolution. This is unsuprising, given that we probe only the first 
lobe of the Fourier transform of the observed source structure, 
which is, therefore, only partially resolved. 

The uncertainties in the amplitude and phase for each spec- 
tral channel were computed from their standard deviations dur- 
ing the averaging of the exposures. Fractional uncertainties of 
only 1 - 2% in the observed visibilities were obtained for all 
channels and baselines. 

2.1. Visibility amplitude and wavelength calibration 

A standard calibration procedure, to assign the corresponding 
wavelength to each spectral channel of the AMBER detector, 
is not currently available. This is a concern, especially for the 
mid-resolution and low-resolution modes, where it is difficult, or 
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UTC 


Seeing 


Air mass 


Coll. time 


EO-GO 


GO-HO 


EO-HO 


PA 




(hh:mm) 


(") 




(ms) 


L(m) 


L(m) 


L(m) 


(deg) 


RS Cap 


07:08-08:06 


0.61/0.94 


1.03/1.12 


1.6/2.4 


14.4/15.5 


28.7/31.0 


43.1 /46.5 


119/113 


18 Cap 


06:47-06:56/08:16-08:26 


0.57/1.19 


1.00/1.11 


1.5/2.7 


14.9 / 16.0 


29.7/31.9 


44.6/47.9 


110/123 



Table 1. Details of our observations. For each star, we give the ranges of values obtained during all the exposures. L(m) is the 
projected baseline length (in meters) and PA is the position angle of the baseline (North through East). Coh. time is the coherence 
time estimated from the atmospheric conditions. 



impossible, to use line profiles either from the observed objects 
or from Earth's atmosphere for the calibration. We performed 
the wavelength (A) calibration of our AMBER observations by 
comparing the positions of the minima of the observed CO band 
heads in RS Cap to their values derived from the models de- 
scribed in Sect. 14.11 The use of different stellar models did not 
affect our A calibration at a level higher than our spectral resolu- 
tion. 

In the amplitude calibration, the CHARM2 catalog of high 
angular-resolution stellar-diameter measurements (Richichi, 
Percheron, & Khristoforova l2003T l indicates that the uniform- 
disk K-band diameter of our calibrator star, 18 Cap, is 5.02 mas. 
This uniform-disk diameter was used to calibrate the visibility 
amplitudes by comparing those of 18 Cap to the model predic- 
tions corresponding to a uniform disk of 5.02 mas diameter. This 
amplitude calibration was performed as described in Annex B of 
the AMBER Data Reduction Software User Manual (Revision 
2.1) of the JMMC0 



2.2. Photometry calibration 

In Fig. \lla), we show the normalized spectrum of RS Cap, ob- 
tained from the average of the photometry channels of the three 
ATs. We note that there is a clear discontinuity in the spectrum 
at ~ 2.2 /im (a similar discontinuity is also seen in the spectrum 
of 18 Cap). This discontinuity is caused by the transmission of 
the K-band fiber of the AMBER instrument. Therefore, it can be 
understood as an instrumental multiplicative effect, which biases 
the photometric and interferometric channels in the same way 
for wavelengths shorter than ~ 2.2 yum. An additive contribution 
to these channels would have affected the contrast of the fringes 
relative to the photometry and, therefore, would have translated 
into a changing amplitude of the visibilities, which is not ob- 
served. Since any multiplicative bias in the data has no effect on 
the visibility amplitudes, the atmospheric opacity and any un- 
modeled gain effect in the detector do not affect the visibilities, 
but only the observed spectra. 

If we wish to use the observed spectrum of RS Cap in our 
modeling, it is necessary to calibrate it with the atmospheric (and 
instrumental) spectral transmission profiles. To perform this cal- 
ibration, we proceeded in the following way: 1) we obtained a 
template spectrum of an MOIII star, which is corrected for atmo- 
spheric effects (star BS4371, observecQby Lan^on et al. 2000); 
2) we divided the template spectrum by the observed spectrum 
of 18 Cap, to obtain the inverse of the transmission profile of the 
atmosphere plus the detector (we call this quantity gain and we 
show it in Fig. [TJb)); and 3) we calibrated the spectrum of RS 
Cap, multiplying it by the gain. The normalized calibrated spec- 
trum is shown in Fig.|2la). By applying this calibration strategy. 




2.15 2.2 2.25 2.3 2.35 2.4 
Wavelength (jum) 



2.45 



^ Tittp://www.inari otti . fr/data_processing_amber .html 

^ This spectrum has a resolution of A/AA ~ 1100, similar to that of 
our observations 



Fig. 1. (a) Normalized spectrum of RS Cap, as observed with 
the ATs; (b) inverse of the atmosphere-i-detector spectral trans- 
mission profile. 



we assume that the airmass is similar for target and calibrator, 
which is indeed the case in our observations (see Table[T]). 

We note that there are some atmospheric telluric lines that 
can be clearly seen in the transmission profile shown in Fig.lTfb); 
those at 2.45 (three Unes), 2.44, 2.42 (three Unes blended), 2.37, 
and 2.32 /urn. The effect of the optical-fiber transmission jump at 
2.2 jum is also seen in the gain, together with a slight increase 
as a function of wavelength redward of ~2.35 fj.m due to wa- 
ter vapor in the Earth's atmosphere. We also note that, although 
the CO bands might have slightly different depths in 18 Cap and 
BS437 1 , no clear features are seen in the regions of the CO bands 
for the gain shown in Fig.lTJb). 



3. Results and discussion 

3. 1 . Continuum angular diameter and effective temperature 

In Fig. |2l we show (a) the normalized spectrum of RS Cap, (b) 
the amplitude visibilities for the three baselines (higher visibility 
corresponds to shorter baseline), and (c) the resulting diameter 
estimates obtained from the fit of a uniformly-bright disk to the 
visibilities at each spectral channel. 
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Fig. 2. (a) Spectrum of RS Cap (normalized); (b) calibrated vis- 
ibility amplitudes; (c) uniform-disc angular diameter at each 
AMBER spectral channel. 



The mean uniform-disk angular diameter fitted to the RS Cap 
visibilities in the continuum (i.e., blueward of 2.29 ;um) is 6 - 
7.75 + 0.07 mas. To derive the effective temperature accurately, 
the size estimated with a limb-darkening model must be used, 
instead of that estimated from a uniform-disc model. We note 
that Teff is a function of the criterion adopted for the radius of the 
star In spherical atmospheres, the actual T^ff used to compare to 
the models should be computed using the diameter at which the 
Rosseland opacity is 1 (i.e., the Rosseland diameter, Oross', see, 
e.g.. Sect. 3.4 of Wittkowski, Aufdenberg, & Kervella llUCT . 

We estimated the Rosseland diameter of RS Cap by fitting 
the visibilities in the continuum part of our spectral coverage 
to the Hankel transform (i.e., Fourier transform of a circularly- 
symmetric source) of the monochromatic intensity profile of 
the MARCS atmospheric model described in Sect. 14.11 We show 
this intensity profile, normalized to a Rosseland radius /?ross = 
0Ross/2 = 1, in Fig. [3] If this profile is given as I(r) (where / is 
the intensity normalized to 1 and r is the radial coordinate nor- 



Radius (Norm.) 

Fig. 3. Monochromatic intensity profile, normalized to Rross = 
1, corresponding to the MRCS model described in Sect. 14. II This 
profile was computed in the continuum part of the spectral cov- 
erage of our observations. 



malized to the Rosseland radius, Rross 
to fit the visibilities is 

V(q)= I(r/RR,MKqr)rdr, 



1), then the model used 



(1) 



where Jq is the first-kind Bessel function of order 0, q is the 
baseline length, and Kisa constant to scale r and q to their corre- 
sponding units (e.g., milliarcseconds and megawavelengths, re- 
spectively). The use of this fitting model, instead of the Hankel 
transform of a uniform disc, allows us to take into account 
limb-darkening and, simultaneously, to estimate directly the 
Rosseland radius (Rross in Eq.[TJ. 

The best-fit Rosseland diameter of RS Cap is 0ross = 7.95 + 
0.07 mas. Our size estimate, together with the bolometric flux of 
(2.1 + 0.2) X 10"'' erg cm"^ s"' reported in Richie hi et al. (lT992l l. 
translates into an effective temperature of Teff = 3160 ± 160 K. 

Lunar occultations of RS Cap, recorded with the TIRGO 
telescope at 2.2 //m, resulted in a uniform-disk diameter estimate 
of 7.75 +0.67 mas (Richichi et al. ll992l l. This diameter, together 
with the bolometric measurements performed by the same au- 
thors, implies an effective temperature of ~ 3200 K, higher than 
the temperature derived from the calibration of Bessel, Castelli 
& Plez ( 1998 ) (~ 2900 K) and that derived from the calibration 
of Scalo (Il976l l (~ 2600 K), but in agreement with the calibra- 
tion given in Ridgway et al. (J980i ). A later re-analysis of the 
Lunar-occultation events (Richichi et al. 1 19991 1. to estimate the 
Rosseland angular diameter, resulted in an even higher effec- 
tive temperature for this star (3481 + 177 K). Dyck, van Belle 
& Thomson ( 119981 1 estimated a size of 7.0 + 0.8 mas from a 
single-baseline interferometric observation in the K band us- 
ing the Infrared Optical Telescope Array (IOTA). This size also 
translates into a higher effective temperature for RS Cap (around 
4000 K). Our estimated diameter of RS Cap is compatible with 
that reported in Richichi et al. (1992) from lunar occultations 
and that reported in Dyck, van Belle & Thomson ( 1998) from 
IOTA observations, although the precision in our estimate is an 
order of magnitude higher. 

3.1 .1 . A word of caution about the use of FINITO 

We note that the use of FINITO in our observations may bias the 
absolute amplitude calibration of the visibilities of each baseline, 
since the percentage of lost fringes during the DIT is not taken 
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Fig. 4. Visibility amplitudes between 2.13 and 2.29 fim as a func- 
tion of baseline length. Black, observations; red, model predic- 
tions with a uniform-disk diameter of 7.75 + 0.07 mas (lower 
line) and 7.75 - 0.07 mas (upper line). 



into account in the process of frame integration of each expo- 
sure. Any possible bias may be different for calibrator and tar- 
get, and translate into a different absolute amplitude calibration 
for each star. As a result, our size estimate of RS Cap might be 
(slightly) biased even after performing the calibration described 
in Sect. 12. II However, the spectral resolution of our observations 
is much higher than that used in Richichi et al. ( 1992 ) and Dyck, 
van Belle & Thomson ( 1998) (who, indeed, only measured the 
visibility at one projected baseline). The higher spectral resolu- 
tion of AMBER allows us to precisely obtain the behavior of vis- 
ibility amplitudes through relatively wide regions of the Fourier 
plane. This sampling of visibility amplitudes through the observ- 
ing band encodes information about the source size, which is in- 
dependent of the absolute calibration of the visibilities. In Fig. 
m we show the visibility amplitudes in the continuum (i.e. blue- 
ward of 2.29 yum) as a function of baseline length. In the figure, 
we also show the model predictions using a uniform-disk diam- 
eter of 7.75 + 0.07 mas. It can be seen that a disk model with this 
diameter satisfactorily fits the visibilities of all baselines through 
the whole spectral coverage in the continuum. This result gives 
us confidence about the amplitude cahbration of the RS Cap vis- 
ibilities. 



3.2. CO first overtone bands: A > 2.3 fim 

In Fig. |2c), there are two clear aspects in the region of CO first 
overtone bands (A > 2.3yL(m) worth noticing. On the one hand, 
we see an increase of ~ 10% in size for all the CO band heads, 
compared to the sizes in the continuum. On the other hand, there 
is a linear trend of increasing size with observing wavelength, 
beginning at 2.29 yum. Owing to this latter effect, the size in the 
continuum at 2.45 /vm is ~ 12% larger than the size in the con- 
tinuum for the wavelengths shorter than 2.29 fim. 

Large size increases redward of 2.29 fim have been observed 
in Mira stars (Woodruff et al. 120091 1. which are interpreted as 
being caused by water-vapor extended zones. However, in those 
cases, the size variations reported are 3-4 times larger than that 
observed in RS Cap. In Sect. 14.21 we discuss how our model of 
water-vapor envelope can explain the observed increasing size 
in RS Cap redward of 2.29 fim. 



4. Modeling 

4.1. MARCS model atmosphere 

We compared the observed size effects in RS Cap with differ- 
ent model atmospheres computed with the MARCS code (e.g. 
Gustafsson et al. 2008 ). The MARCS code computes the hydro- 
static atmospheric structure in radiative and convective equilib- 
rium for spherical atmospheres with molecular and atomic lines 
taken into account using the opacity sampling method. Each 
model in spherical geometry (suitable for AGB stars) is speci- 
fied by the effective temperature (Tetf), surface gravity (logg), 
stellar mass (M), microturbulent velocity, and chemical compo- 
sition. 

For Teff , we adopted 3200 K, the nearest possible value in 
MARCS models to our estimated 3160K. We also assume the 
"moderately CN-cycled" chemical composition (C/N =1.5 and 
'^C/'^C = 20), which is usually found in stars after the red giant 
branch, and [Fe/H] - 0.0. The stellar mass can be estimated by 
comparing with theoretical stellar evolutionary tracks, once T^^ff 
and the luminosity are known. While the parallax of RS Cap is 
not precisely measured (the relative error is on the order of 70%, 
so the error in the luminosity can be as large as an order of mag- 
nitude). Winters et al. (120031 1 estimated the luminosity of RS Cap 
to be log L/Lq = 3.82, by applying the period-luminosity re- 
lation for semi-regular variables given by Feast (19960. From 
Fig. 5 of Feast (1996 ), we estimated an uncertainty of +0.25 in 
the bolometric absolute magnitude of RS Cap, which translates 
into an error of up to 20-25% in the luminosity, much smaller 
than that related to the Hipparcos parallax. The luminosity from 
Winters et al. ( 120031 1. together with the above Tetj, places RS 
Cap close to the evolutionary track of a 2 Mq star with Z = 0.02, 
computed by F. Herwi^. Therefore, we adopted 2Mo for RS 
Cap. 

We note that the luminosity estimated by Winters et al. 
(2003) for RS Cap is much higher than that corresponding to an 
RGB star with 2 Mq (which is only log L/Lq ~ 3). Therefore, ac- 
cording to its luminosity, RS Cap would be clasified as an AGB 
star, as it is indeed assumed in the P-L relation by Feast ( .1996b . 
The luminosity estimated in this way for RS Cap is, thus, self- 
consistent in the framework of the P-L relation of semi-regular 
AGB stars. 

From the luminosity and T^ff, we estimated a stellar ra- 
dius of ~ 260 Rq, which, combined with our angular-size esti- 
mate, translates into a distance of about 310pc. Hence, from the 
adopted stellar mass of 2 M© and the estimated radius, we de- 
rive log^ ~ -0.1. We adopted logg = 0.0 for the selection of 
the MARCS models. For the microturbulent velocity, we used the 
MARCS models with 2 and 5 kms"'. However, the effects of this 
parameter on the visibility are very small. 

We downloaded the models with the above parameters from 
the MARCS websiteQ Using the temperature and pressure strati- 
fications of the model, we computed the monochromatic inten- 
sity profiles and the synthetic spectrum for the spectral range 
of our AMBER observations with a fine wavelength mesh of 
1.4km s ', sampling several points over each line profile. The 
CO ('^C'^O, '^C'^O, i2ci**0, and '^C' ^O) an d HjO lines were 
taken from the line lists of Goorvitch ( 119941 ) and Partridge & 
Schwenke ( 119971 ), respectively. The monochromatic visibilities 
were computed from the Hankel transform (i.e., Fourier trans- 

Although the period-luminosity relation is a matter of debate, as 
Winters et al. ( 2003 ) mention. 
' http : // astrowww . phy s . uvic . ca/ ~£herwig"| 
|http : //marcs . astro . uu . se | 
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form of a circularly-symmetric distribution) of the monochro- 
matic intensity profiles (computed as described in Ohnaka et 
al. 2006) and then convolved to the spectral resolution of our 
AMBER measurements. 

In Quirrenbach et al. (THOVt, limb-darkened photospheric 
models were used to generate equivalent disk sizes. These sizes 
were then compared to the disk sizes fitted to the Mklll visibili- 
ties. In our case, this approach would have had the disadvantage 
of being sensitive to scaling factors between different baselines 
(i.e., to the absolute amplitude calibration) and would have de- 
pended on the use of an intermediate model (a uniform disk) to 
compare observations to models. To avoid those drawbacks, we 
chose a new approach to the modeling of the RS Cap visibili- 
ties, which minimizes the effect of possible calibration artifacts 
between baselines and, in addition, avoids the use of an inter- 
mediate fit to a disk. The synthetic visibilities obtained from the 
MARCS models were computed at the points of the Fourier plane 
corresponding to each exposure in our observations. We then 
averaged the synthetic visibilities over all the exposures, in the 
same way as we did with the real data. In this way, we directly 
compare the observed visibilities to the model predictions. 

We show in Fig. |6]a comparison between observed visibil- 
ities (black lines) and those computed from the MARCS model 
using Teff = 3200 K, M = 2Mo, and log^ = (blue Hnes). 
The match between observations and model predictions in the 
continuum (i.e., for wavelengths shorter than ~ 2.3 fim) is very 
good. 

The reduced for this model isx^ed ~ 18.2. The relatively 
large depth of the features of the CO band heads cannot be re- 
produced by the model. The discrepancy remains even if we 
lower the mass to values as low as M = 1 (models with 
higher masses predict more compact atmospheres and, there- 
fore, smaller size effects in the CO band heads). Additionally, 
the observed visibility amplitudes at wavelengths longer than 
~2.29//m are systematically lower than the model predictions; 
this difference between model and observations is larger at 
longer wavelengths. In other words, the effective size of the star 
at wavelengths longer than ~ 2.29 jum is systematically larger 
than the effective size derived from the MARCS model. 

4.2. An additional water-vapor envelope 

An extra contribution must be added to the MARCS models to fit 
the visibilities in the observed spectral range. We can reproduce 
the lower visibility amplitudes (i.e., the larger effective sizes) be- 
yond 2.29 yum by adding a contribution to the stellar opacity due 
to water vapor around the star The opacity corresponding to the 
wide absorption of water vapor centered at -2.7 ^va may account 
for the larger angular sizes observed at the longer wavelengths. 

We added a simple water- vapor model envelope to the MARCS 
pressure/density profiles, consisting of a narrow spherical shell 
with a radius of 2 and a width of 0. 1 (both in units of the stel- 
lar radius), a column density of 10^' cm"^, and a temperature of 
1500 K. Using these parameters, we are able to fit the amplitude 
visibilities at /I > 2.29 yum, as we show in Fig|6](red lines). With 
this model, we obtain = 8.1. This large reduced is par- 
tially due to the discrepancy in the CO bands, which is discussed 
in Sect. 14.31 Our model follows the general trend of decreas- 
ing visibility amplitudes for wavelengths longer than 2.3 /^m, al- 
though underestimates the visibility amplitudes at wavelengths 
longer than 2.4yum. Changing the column density by a factor of 
ten results in an increase x^g^ by ^ factor 2 - 3; a slight mod- 
ification of the water-vapor temperature (~ 15%) increases 




2.15 2.2 2.25 2.3 2.35 2.4 2.45 
Wavelength (jum) 



Fig. 5. Calibrated spectrum of RS Cap (black), synthetic MARCS 
spectrum (blue), and synthetic spectrum with the additional 
water- vapor envelope (red). 



by a factor ~ 2; and changing the size of the shell by a factor of 
two increases the x^^^ by a factor ~ 2. 

Although the water-vapor envelope enhances the contrast of 
the CO features in the visibilities, thus improving the fit also in 
the CO band heads, the observed depths in the band heads are 
still larger than those of the model for some band heads (~ 50% 
for the band head at 2.29 /im). Furthermore, increasing (decreas- 
ing) the stellar mass results in smaller (larger) depths of the visi- 
bility amplitudes in the CO band heads, but decreasing the mass 
to a value as low as 1 M© does not translate into CO visibility 
amplitudes comparable to those of the observations, even after 
the addition of the water envelope (the x^^^ using 1 Mq is only 
7.2% lower than that using 2 Mq). 

In Fig. |5j we show the calibrated spectrum of RS Cap 
(black) and, superimposed, the synthetic spectra obtained from 
the MARCS model only (blue) and the MARCS model with an addi- 
tional water-vapor envelope (red). Both model spectra are very 
similar, and only differ for the longest wavelengths, where the 
model with water envelope predicts a slightly lower flux. Since 
the fit of both models to the observed spectrum is of similar qual- 
ity, it is difficult to extract any additional conclusion about the 
water-vapor envelope from the spectrum alone. 

4.2.1. Water-vapor envelopes in AGBs. The problem of RS 
Cap 

The presence of dense water vapor envelopes in semi-regular or 
irregular (i.e., non-Mira-type) AGB stars was first revealed by 
Tsuji et al. (119971) . who unmasked the 2.1 water vapor band 
originating in the dense molecular layers extending to ~2 stel- 
lar radii. On the basis of infrared interferometric observations, 
Mennesson et al. ( I2002i) found an increase in the angular size of 
30-60% between the K and L' band in five semi-regular AGB 
stars, which the authors interpreted as being due to the emis- 
sion from the extended molecular layers. However, they used fil- 
ters which are incapable of spectrally resolving the water vapor 
bands. Our AMBER measurements are the first study to spatially 
and spectrally resolve the water vapor emission from the dense 
molecular layers for a semi-regular AGB star. 

We note that a similar increase in the angular diameter red- 
ward of 2.3jL/m is also found in Mira stars. Based on /T-band 
interferometry using narrow-band filters with a spectral resolu- 
tion of ~20, Perrin et al. ( I2004I I showed that the angular diam- 
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eters of Mira stars are larger at 2.0 and 2.4 yum than at 2.2 yum. 
Wittkowski et al. (2008 ) observed the Mira star S Ori from 1.29 
to 2.32yL(m using AMBER with a spectral resolution of 35 and 
found that the angular size increases redward of 2.2 jum. The 
wavelength dependence of the angular size of Mira stars from 
the near- to mid-infrared is primarily governed by the opacity of 
water vapor, which is abundant in the atmosphere of Mira stars 
(e.g., Ohnak a 2004, Wittkowski et al. 2008). The model by Tej 
et al. ( I2003I I shows that a shell with an enhanced water vapor 
abundance forms behind a shock front propagating outward in 
the extended atmosphere of Mira stars. Comparison between in- 
frared interferometric observations of Mira stars and dynamical 
models lends support to this picture (e.g., Ohnaka et al. 120061 
Wittkowski et al. 120071 ). 

The shock wave in Miras is generated by the periodic, large- 
amplitude pulsation. Therefore, the shock in non-Mira stars 
with much smaller variability amplitude is expected to be much 
weaker. However, despite the striking difference in the variabil- 
ity amplitude (AV - 6-9 and 0.5-2 for Miras and RS Cap, re- 
spectively), the semi -regular AGB star RS Cap produces a warm 
water vapor envelope similar to that of Miras. The column den- 
sity, temperature, and radius of the water vapor envelope derived 
for RS Cap are similar to those derived for Mira stars (Ohnaka 
120041 1 ■ Not only the properties of the water-vapor envelope, but 
also the mass-loss rate of RS Cap (~ 1 x 10"^ M© yr"') is compa- 
rable to that of Miras (e.g., 5x10"'' Mq yr"' for S Ori) (estimates 
from Winters et al. 2003). 

The dynamical models with small-amplitude pulsations, re- 
ported in Winters et al. (120031 ). produce outflows slower than 
5 km s"' with mass-loss rates lower than 3 x 10"^ Mq yr"'. 
However, the expansion velocity measured for RS Cap is ~ 1 km 
s"' (Winters et al. 2003 ). Therefore, it is not clear whether this 
kind of small-amplitude pulsation model applies to RS Cap. At 
the moment, there is no satisfactory explanation for the origin of 
the warm water envelope in RS Cap. 



4.3. An additional CO envelope 

It can be seen in Fig. |6]that the observed visibility amplitudes at 
the CO band heads are much lower than the model predictions, 
either with or without a water-vapor envelope. This discrepancy 
is especially large for the shortest baseline (i.e. £0 - GO) at the 
2.29/im band head. The low visibility amplitudes at the CO band 
heads provide evidence of an extended envelope of CO around 
the star For instance, if ~20% of the emission in the CO band 
head at 2.29 fim were to originate from a very extended envelope 
of, say, ~20 mas radius or larger, it would map into a decrease in 
the visibility amplitude that would roughly match the observa- 
tions in that part of the spectrum. Unfortunately, we do not have 
enough data to be able to characterize well this possible extra 
component of CO absorption. 

5. Conclusions 

We have observed the AGB star RS Cap with VLTI/AMBER 
in the K band with medium-resolution mode (506 channels be- 
tween 2. 13 and 2.47 fim). We have estimated a Rosseland diam- 
eter of 7.95 + 0.07 mas in the continuum, which translates into 
an eff'ective temperature of 3 1 60+1 60 K. The apparent size of 
the star increases monotonically by ~ 12% between 2.29 and 
2.47 /im. We have detected lower than expected visibility am- 
plitudes in all the CO band heads observed. These lower ampli- 
tudes translate into larger apparent sizes of the star in the CO 
band heads. 
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Fig. 6. Black: AMBER visibihties of RS Cap. Blue: model 
visibilities computed using a MARCS pressure-density profile 
(see text). Red: model visibilities computed using the MARCS 
pressure-density profile and an additional water-vapor envelope 
(see text). 



Using the pressure/density model profiles obtained with the 
MARCS code, we have been able to generate synthetic visibilities 
and compare them directly to our AMBER observations. The fit 
is rather good, although we are unable to reproduce the low vis- 
ibility amplitudes in the CO band heads using a mass for RS 
Cap estimated from the theoretical models of stellar evolution. 
In this sense, our situation resembles that of Quirrenbach et al. 
(120011) . who observed interferometrically a sample of cool-giant 
stars in the TiO absorption band at 7 12 nm. The discrepancy be- 
tween models and observations in the CO band heads might be 
resolved either by using a much lower mass for RS Cap (below 
1 Mq, thus in contradiction with the stellar-evolution models) or 
by additional unmodeled effects in the base of the stellar wind (a 
transition zone with prominent absorption in the CO bands). 

To fit the observations at wavelengths longer than 2.29yL(m, 
we have found that an ad hoc narrow spherical water-vapor en- 



8 



I. Marti- Vidal et al: AMBER observations of RS Cap 



velope around the star, similar to the models used in Perrin et al. 
(I2004I I. must be added. We have modeled this envelope with a 
temperature of 1500K, a size twice that of the star, a width of 
0.1 times the stellar radius, and a column density of 10^' cm"^. 

Finally, there is a hint of an extended CO envelope around 
the star, based on the low visibility amplitudes at the CO band 
heads in the shortest baseline. 
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